Chaotic exit to inflation: the dynamics of pre-inflationary universes 
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We show that anisotropic Bianchi type-IX models, with matter and cosmological constant 
have chaotic dynamics, connected to the presence of a saddle-center in phase space. The 
topology of cylinders emanating from unstable periodic orbits about the saddle-center 
provides an invariant characterization of chaos in the models. The model can be thought 
to describe the early stages of inflation, the way out to inflation being chaotic. 



Some important features in inflationary dynamics arise whenever anisotropy 
is present even in the form of small perturbations. This point has not been em- 
phazised yet in the literature of inflation. The conjunction of the cosmological 
constant, anisotropy and matter fields implies the existence of a saddle-center E in 
phase spacil. As consequence we have a complex dynamics, based on structures as 
homoclinic cylinders, which emanate from unstable periodic orbits in the neighbor- 
hood of E. These cylinders will cross each other transversally producing a chaotic 
dynamics, analogous to the breaking of homoclinic curves in Poincare homoclinic 
phenomeuEU. This structure constitutes an invariant characterization of chaos. 

We consider anisotropic Bianchi IX cosmological models characterized by two 
scale factors A{t) and B{t), with dust plus a cosmological constant A. The dynamics 
of the models is governed by the Hamiltonian 

Pa Pa A p2 43 
H(A,B,Pa,Pb) = 4- + 2 A ^ ~2AAB'^ - Eo = 0, (1) 

where Pa and Pb are the momenta canonically conjugated to A and B, respectively, 
and Eq is a constant proportional to the total energy of the models. The associated 
dynamical system has one saddle-center E in the finite region of the phase space 
with associated energy Eq = Ecrit ~ ^-nd eigenvalues Ai,2 = ± — , A3 4 = 
±2 iEcrit ■ The phase space has two critical points at infinity, corresponding to 
the de Sitter solution one of them acting as an attractor. The dynamical system 
generated by (1) admits the invariant manifold ^A {A — B, Pa = Its phase 

portrait is depicted in Fig. 2, where the curves represent isotropic universes. The 
critical point E is contained in A4, and the separatrices S constitute boundaries 
between isotropic models that collapse and escape to the de Sitter configuration. 

Due to its saddle-center character, the motion in a small neighborhood of the 
critical point E is separable into hyperbolic and rotational motions, with respective 



energies approximately conserved. The rotational motion corresponds to periodic 
orbits in a linear neighborhood of E (cf. Fig 1(a)). For zero energy of the hyperbolic 
motion, we have the linear stable Vg and unstable Vu one-dimensional manifolds of 
Fig. 1(b). The separatrices S in the invariant manifold Ai are actually the non- 
linear extension of Vg and Vu ■ The direct product of a periodic orbit with Vg and 
Vu generates in the linear neighborhood of E the structure of stable and unstable 
cylinders, which coalesce into the periodic orbit for large positive and negative times, 
respectively (Fig. 1(c)). A general orbit which visits the neighborhood of E has 
an oscillatory approach to the cylinders, the closer as the energy of the hyperbolic 
motion goes to zero. The outcome of this oscillatory regime will be to collapse or to 
escape to the de Sitter attractor depending on the energy of the hyperbolic motion. 

Away from a linear neighborhood of E^ the linear approximation is no longer 
valid. Higher order terms become important for the dynamics, and the non-integra- 
bility of the system results in the distortion and twisting of the cylinders. The stable 
cylijader and the unstable one will cross each other transversally, producing chaotic 
setffl in the phase space; orbits with initial conditions taken on these chaotic sets 
will be highly sensitive to small perturbations on the initial conditions. 

The phase space under consideration is not compact, and we will actually iden- 
tify chaotic behaviour associated to the possible asymptotic {t oo) outcomes of 
orbits, namely, escape to the de Sitter attractor at infinity or collapse after a burst of 
initial expansion. We assume A = 0.25, so that the critical point E is characterized 
hy A = B = 1.0, Pa Pb = and Ecrit = 1.0. All calculations were made using 
the package PoincarQ. Let Sq be a point belonging to the separatrix {Eq = 1.0), 
with coordinates A = B = 0.4, Pb ^2Pa^ 1.3576450198. Around this point, we 
construct a 4-dim sphere in the phase space with arbitrary small radius Rq, for in- 
stance Rq = 10^'"'. The values of {A, B, Pa, Pb) are taken in energy surfaces which 
have a non-empty intersection with this sphere. Such energy surfaces are those for 
which the range of energy, Eq, about Eq = Ecrit = 1-0 is of the order of, or smaller 
than the radius Rq. Physically we are considering expanding cosmological models 
with small anisotropic perturbations in a pre-inflationary phase. 

After several experiments with the above set of initial conditions, we note that, 
as expected, two possible outcomes arise: collapse, or expansion into the de Sitter 
configuration. The main result of this paper is to show that, for a determined 
interval of energy 5 E*, this outcome is chaotic. Indeed, according to our numerical 
work, we find out that, for each sphere of initial conditions, there always exists a 
non-null interval 5 E* for which orbits either collapse or escape, in an indeterminate 
outcome. In Fig. 3 this behaviour is showed for spheres of initial conditions with 
radius i?o = 10^^. An empirical relation between the gap 5 E* and the radius i?o is 
obtained, 5 E* oc Rq. In other words, any infinitesimal fluctuation of a given initial 
condition in the sphere (for energies in the interval S E* can lead to an indeterminate 
outcome, that is, collapse or escape. This is evidence of chaos, as a consequence 
of the crossing of stable and unstable cylinders, emanating from unstable periodic 
orbits about E. This topological structure is actually an invariant characterization 
of chaos. Finally we remark that the above behaviour is not restricted to initial 
conditions taken in small neighborhoods of points of the separatrix. Generally, 
any sets of initial conditions taken in an arbitrary neighborhood of the invariant 
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manifold A4, which result in orbits that visit a small neighborhood of E, display 
the above chaotic behaviour. 




Figure 1: (a)Periodic orbits of the Hamiltonian in teh linear approximation. (b)The linear unstable 
Vu and stable Vs one-dimensional manifolds. (c)Stable and unstable cylinders manifolds emanating 

from the periodic orbit r. 




Figure 2: Integral curves on the invariant man- Figure 3; Chaotic exit to inflation: outcome of 
ifold A = B, Pa = PbI"^ for 7=1 (dust). 50 orbits choosen in a sphere of -R = 10~® about 

a point 5o on the separatrix. 
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